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BESCHREIBIJNfi- 

Bildgebendes Verfahren mit gefilterter Ruckprojektion 

Die Erfindung betrifft ein bildgebendes Verfahren, bei dem ein Projektionsdatensatz 
eines zu rekonstruierenden Untersuchungsbereichs dutch Akquirieren von Projektionen 
aus unterschiedlichen Projektionsrichtungen, insbesondere mit Hilfe eines Computer- 
tomographen, erzeugt und aus diesem Projektionsdatensatz ein Bild des Untersuchungs- 
bereichs mit Hilfe einer gefilterten Ruckprojektion rekonstruiert wird. Die Erfindung 
bezieht sich auBerdem auf einen Computertomographen zur Durchfuhrung des Ver- 
fahrens sowie auf ein Computerprogramm zur Steuerung des Computertomographen. 

In bekannten Verfahren der eingangs genannten Art konnen Projektionsdatensatze 
bspw. mit einem Positronen-Emissions-Tomographen (PET), einem Einzel-Photon- 
Emissions-Computertomographen (Single Photon Emission Computed Tomography - 
SPECT), einem C-Arm-System oder einem Computertomographen erzeugt werden. 

Bei einem Computertomographen wird der Untersuchungsbereich von einem Strahlen- 
biindel durchstrahlt, von der mtensitat des Strahlenbundels jenseits des Untersuchungs- 
bereichs abhangige Messwerte werden detektiert und ein Computertomographie-Bild 
(CT-Bild) wird aus diesen Messwerten mit Hilfe einer gefilterten Ruckprojektion 
rekonstruiert. Die raumliche Verteilung der Absorption bzw. der Schwachung der 
Strahlung in dem Untersuchungsbereich kann mittels einer gefilterten Ruckprojektion 
rekonstruiert werden. 

Die gefilterte Ruckprojektion hat den Vorteil, dass sie einen im Vergleich zu anderen 
Rekonstruktionsverfahren geringen Rechenaufwand erfordert, zu Bildern mit sehr guter 
Qualitat fuhrt, eine Rekonstruktion beliebiger Unterbereiche des Untersuchungsbereichs 
und ein gleichzeMges Verarbeiten unterschiedlicher Projektionen bzw. Messwerte, d.h. 
eine parallele Verarbeitung von Messwerten, ermSglicht. 
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Problematisch ist bei der gefilterten Ruckprojektion das Ermitteln desjenigen Filters, 
der fiir bestimmte Projektionsgeometrien (beim Computertomographen bspw. definiert 
durch den Verlauf der Strahlen einer Projektion) zu Bildem mit optimaler Bildqualitat 
fUhrt Der entsprechende Filter kann bspw. durch Versuche mit unterscbiedlichen 

5 Filtern oder mit sehr grofiem Aufwand analytisch ermittelt werden. Da das analytische 
Bestimmen des Filters sehr schwierig ist, sind fur viele Projektionsgeometrien (beim 
Computertomographen bspw. fiir Strahlenverlaufe, die durch helixformige Trajektorien 
mit beliebigem Tischvorschub pro Rotation und kegelformigem Strahlenbundel erzeugt 
werden) FUter bestimmt worden, die die jeweilige Projektionsgeometrie nur unzu- 

10 reichend und naherungsweise beriicksichtigen, was zu artefaktreichen rekonstruierten 
Bildern fuhrt 

Es ist daher die Aufgabe der vorhegenden Erfindung ein bildgebendes Verfahren, einen 
Computertomographen und ein Computerprogramm anzugeben, bei denen zur Rekon- 
15 straktion eines Bildes ein fur die jeweilige Projektionsgeometrie optimaler Filter 
ermittelt und verwendet wird, so dass fur jede Projektionsgeometrie Bilder von sehr 
guter Bildqualitat erzeugt werden. 

Diese Aufgabe wird erfindungsgemaB gelost durch ein bildgebendes Verfahren mit den 
20 Schritten: 

a) Erzeugen eines Projektionsdatensatzes eines zu rekonstruierenden Unter- 
suchungsbereichs durch Akquirieren von Projektionen aus unterscbiedlichen 
Projektionsrichtungen, wobei jeder Projektion eine Projektionsrichtung und 

25 jeder Projektionsrichtung eine Projektionsgeometrie zugeordnet ist, 

b) Bereitstellen eines Iterationsoperators, der abhangig ist von einer ersten Projek- 
tionsrichtung, einer zweiten Projektionsrichtung und einem der zweiten Projek- 
tionsrichtung zugeordneten Filteroperator, und Ermitteln jeweils eines Filter- 

30 operators fur jede Projektionsrichtung, wobei die Filteroperatoren so angepasst 
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werden, dass eine Mittelung fiber alle Iterationsoperatoren, die dieselbe erste 
Projektionsrichtung aber unterschiedliche zweite Projektionsrichtungen 
aufweisen, im wesentlichen einen Identitatsoperator ergibt, und wobei bei 
Anwendung eines Iterationsoperators auf eine Projektion folgende Operationen 
5 ausgefuhrt werden: 

- Rfickprojizieren der Projektion in einer Projektionsgeometrie, die der ersten 
Projektionsrichtung zugeordnet ist, 

1 0 - Vorwartsprojizieren der ruckprojizierten Werte in einer Projektions- 

geometrie, die der zweiten Projektionsrichtung zugeordnet ist, 

- Filtern der vorwarts projizierten Werte mit einem Filteroperator der zweiten 
Projektionsrichtung, 

15 

- Rfickprojizieren der gefilterten Werte in einer Projektionsgeometrie, die der 
zweiten Projektionsrichtung zugeordnet ist, 

- Vorwartsprojizieren der gefilterten, ruckprojizierten Werte in einer 
Projektionsgeometrie, die der ersten Projektionsrichtung zugeordnet ist, 

c) Rekonstruieren eines Bildes des Untersuchungsbereichs aus den Projektionen 
mit Hilfe einer gefilterten Rfickprojektion, wobei auf jede Projektion der 
Filteroperator der entsprechenden Projektionsrichtung angewendet wird. 

Jede Projektion weist mehrere Messwerte auf, die der jeweiligen Projektionsrichtung 
zugeordnet werden konnen. Die jeder Projektion zugeordnete Projektionsgeometrie 
beschreibt, wie eine Vorwarts- und Rfickprojektion fur diese Projektion bzw. die 

entsprechende Projektionsrichtung ausgefuhrt werden kann. Dies wird unten 
beispielhaft fur einen Computertomographen nahcr erlautert 
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Das Riickprojizieren kann durch einen Ruckprojektions-Operator beschrieben warden, 
der eine Projektion von einem Projektionsraum in den Untersuchungsbereich projiziert. 
Ein Projektionsraum kann bspw. bei einem Computertomographen ein Raum sein, der 
durch die Detektorflache und die Strablenquellenposition aufgespannt wird. Das 

5 Vorwartsprojizieren kann durch einen Vomartsprojektions-Operator beschrieben 

werden, der von dem Untersuchungsbereich in den Projektionsraum vorwarts projiziert. 
Alternativ kann auch ein Operator definiert werden, der nacheinander eine Ruck- und 
eine Vorwartsprojektion ausfuhrt. Dieser Operator bildet eine Projektion von dem 
Projektionsraum in denselben ab und kann die entsprechende Kombination aus Ruck- 

10 und Vorwartsprojektions-Operatoren ersetzen. 

Die im Schritt b) beschriebene Ermittlung eines Filteroperators fuhrt zu einem an die 
jeweilige Projektionsgeometrie angepassten, optimalen Filter, der nach Rekonstruktion 
mittels einer gefilterten Riickprojektion zu Bildem mit sehr guter Qualitat tuhrt. Dies 
15 fuhrt insbesondere bei Projektionsgeometrien, fur die bisher nur Naherungslosungen fur 
Filter bekannt shid, zu einer im Vergleich zum Stand der Technik verbesserten 
BildquaUtat 

Die in Anspruch 2 beschriebene iterative Ausfuhrungsform ermoghcht die Bestimmung 
20 optimaler Filteroperatoren aufeine mathematisch relativ einfache Weise, wodurch der 
Rechenaufwand weiter reduziert wird. 

In Anspruch 3 wird eine Ausfuhrungsform beschrieben, in der ein CT-Bild aus einem 
Projektionsdatensatz, der mit Hilfe eines Computertomographen erzeugt wird, rekon- 
25 struiert wird, das eine im Vergleich zu bekannten Computertomographie-Verfahren 
verbesserte Bildqualitat aufweist. 

In Anspruch 4 wird eine Ausgestaltung beschrieben, in der der Filteroperator fur alle 
Projektionsrichtungen gleich ist. Dies fuhrt zu einem reduzierten Rechenaufwand bei 
30 unverandert guter Bildqualitat. 
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Die erfindungsgemaBe Rekonstruktion des Untersuchungsbereichs aus einem 
Projektionsdatensatz, der nach Anspruch 5 mit einem PET oder einem SPECT erzeugt 
worden ist, fiihrt zu Bildern von sehr guter Qualitat. 



5 Die Diskretisierung zumindest einiger Operatoren nach Anspruch 6 ermoglicht eine 
unkomplizierte, relativ einfache Implementierung des erfindungsgemaBen Verfahrens. 

Die Darstellung zumindest einiger der Operatoren durch lineare Funktionale nach 
Anspruch 7 fiihrt zu einer exakten Berechnung der Filteroperatoren und damit zu einer 
1 0 weiter verbesserten Bildqualitat. 



Ein Computertomograph zur Durchfuhrung des erfindungsgemaBen Verfahrens ist in 
Anspruch 8 beschrieben. 

15 Anspruch 9 definiert ein Computerprogramm zur Steuerung eines Computertomo- 
graphen nach Anspruch 8. 

Die Erfindung wird im Folgenden anhand der Zeichnungen naher erlautert. Es zeigen: 

20 Fig. 1 einen Computertomographen, mit dem das erfindungsgemaBe Verfehren 
ausfuhrbar ist, 

Fig. 2 ein Ablaufdiagramm des erfindungsgemaBen Verfehren und 
25 Fig. 3 ein Ablaufdiagramm der erfindungsgemaBen Bestimmung von Filteroperatoren. 



Der in Fig. 1 dargestellte Computertomograph umfasst eine Gantry 1, die urn eine 
parallel zur z-Richtung des in Fig. 1 dargestelltenKoordinatensystems verlaufende 
Rotationsachse 14 rotieren kann. Dazu wird die Gantry 1 von einem Motor 2 mit einer 
30 vorzugsweise konstanten, aber einstellbaren Winkelgeschwindigkeit angetrieben. An 
der Gantry 1 ist eine Strahlenquelle S befestigt, beispielsweise ein Rontgenstrahler. 
Dieser ist mit einer Kollimatoranordnung 3 versehen, die aus der von der Strahlenquelle 
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S erzeugten Strahlung ein kegelfonniges Strahlenbiindel 4 ausblendet, dh. ein 
Strahlenbiindel, das sowohl in z-Richtung als auch in einer dazu senkrechten Richtung 
(dh. in einer zur Rotationsachse senkrechten Ebene) eine von Null verschiedene, 
endliche Ausdehnung hat. Alternativ konnte auch ein Scherffirmiger Strahl verwendet 
werden. 

Das Strahlenbiindel 4 durchdringt einen zylinderformigen Untersuchungsbereich 13, in 
dem sich ein Objekt (nicht dargestellt) befindet. 

Nach dem Durchsetzen des Untersuchungsbereichs 13 trifft das Strahlenbiindel 4 auf 
eine an der Gantry 1 befestigte Detektoreinheit 1 6 mit einer Detektorflache, die eine 
Vielzahl von Detektorelementen umfosst, die in dieser Ausfuhrungsform in Zeilen und 
Spalten matrixformig angeordnet sind Die Detektorspalten verlaufen vorzugsweise 
parallel zur Rotationsachse 14. Die Detektorzeilen befinden sich in zur Rotationsachse 
senkrechten Ebenen, in dieser Ausfiihrungsform auf einem Kreisbogen urn die Strahlen- 
quelle S (fokus-zentrierter Detektor). In anderen Ausfuhrungsformen konnen sie aber 
auch anders geformt sein, z. B. einen Kreisbogen urn die Rotationsachse 14 beschreiben 
oder geradlinig sein. Jedes von dem Strahlenbiindel 4 getroffene Detektorelement liefert 
in jeder Position der Strahlenquelle einen Messwert fur einen Strahl aus dem Strahlen- 
biindel 4. Wenn in anderen Ausfuhrungsformen ein facherformiges Strahlenbiindel 
verwendet wird, so konnte die Detektoreinheit auch nur eine Detektorzeile auftveisen. 

Der mit a mgx bezeichnete Of&iungswinkel des Strahlenbiindels 4 bestimmt den 
Durchmesser des Objektzylinders, innerhalb dessen sich das zu untersuchende Objekt 
bei der Akquisition der Messwerte befindet Dabei ist der Of&ungswinkel als der 
Winkel definiert, den ein Strahl, der in einer zur Rotationsachse 14 senkrechten Ebene 
am Rande des Strahlenbiindels 4 liegt, mit einer durch die Strahlenquelle S und die 
Rotationsachse 14 definierten Ebene einschlieBt Der Untersuchungsbereich 13 bzw. das 
Objekt oder der Patientenlagerungstisch konnen miltels eines Motors 5 parallel zur 
Rotationsachse 14 bzw. zur z-Achse verschoben werden. Dazu aquivalent konnte aber 
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auch die Gantry in diese Richtung verschoben werden. Wenn es sich um ein technisches 
Objekt handelt und nicht um einen Patienten, kann das Objekt bei einer Untersuchung 
gedreht werden, wahrend die Strahlenquelle S und die Detektoreinheit 16 still stehen. 

5 Mit Hilfe der Motoren 2 und 5 konnen die StrahlenqueUe S und die Detektoreinheit 1 6 
eine Trajektorie relativ zum Untersuchungsbereich 13 beschreiben, die auf einer 
imaginaren Zylinderoberflache verlauft. Diese Trajektorie kann bspw. helixformig 
verlaufen, wenn beide Motoren im Betrieb sind. Wenn dagegen der Motor 5 fur den 
Vorschub in Richtung der Rotationsachse 14 stillsteht und der Motor 2 die Gantry 
10 rotieren lasst, ergibt sich eine kreisformige Trajektorie fur die Strahlenquelle S und die 
Detektoreinheit 16 relativ zum Untersuchungsbereich 13. 

Die von der Detektoreinheit 16 akquirierten Messwerte werden einer Rekonstruktions- 
einheit, insbesondere einem Rekonstruktionscomputer 10 zugefuhrt, der mit der 

15 Detektoreinheit 16 z. B. iiber eine drahtlose Datenubertragung (nicht dargestellt) 
vefbunden ist. Die Rekoristruktionseinheit 10 rekonstruiert eine Verteilung von 
Schwachungswerten, bspw. eine Absorptionsverteilung im Untersuchungsbereich 13 
und gibt sie, z. B. auf einem Monitor 1 1, wieder. Die beiden Motoren 2 und 5, die 
Rekonstruktionseinheit 10, die Strahlenquelle S und der Transfer der Messwerte von der 

20 Detektoreinheit 16 zur Rekonstruktionseinheit 10 werden von der Steuereinheit 7 
gesteuert. 

In anderen Ausfuhrungsformen konnen die akquirierten Messwerte zur Rekonstruktion 
zunachst einem oder mehreren Rekonstruktionscomputern zugefuhrt werden, die die 
25 rekonstruierten Daten zB. iiber ein Glasfaserkabel an einen Bildverarbeitungscomputer 
weiterleiten. 

Fig. 2 zeigt den Ablauf eines Mess- und Rekonstruktionsverfahrens, das mit dem 
Computertomographen nach Fig. 1 durchgefuhrt werden kann. 

30 
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Nach der tatialisierung im Schritt 101 rotiert die Gantry mit einer Wirdcelgeschwindig- 
keit, die in diesem Ausfuhrungsbeispiel konstant ist. Sie kann aber auch variieren, z.B. 
in Abhangigkeit von der Zeit oder von der StraMenquellenposition. 

Im Schritt 102 wird der Untersuchungsbereich 13 bzw. das Objekt oder der Patienten- 
lagerungstisch parallel zur Rotationsachse 14 verschoben und die Strahlung der 
Strahlenquelle S wird eingeschaltet, so dass die Detektorcinheit 16 die Strahlung aus 
einer Vielzahl von Winkelpositionen ertassen kann. 

Im Schritt 103 werden die Messwerte gefiltert, indem ein Filteroperator auf die 
Messwerte angewandt wird. Das Filtern der Messwerte kann durch folgende Gleichung 
beschrieben werden: 

p t =F iPi . (1) 

Hierbei ist Pi eine der im Schritt 102 akquirierten Projektionen, wobei eine Projektion 
die Menge aller Messwerte ist, die an einer bestimmten Strahlenquellenposition 
akquiriert worden sind. Jeder Projektion kann eine Projektionsrichtung zugeordnet 
werden, die durch die jeweilige Strahlenquellenposition und die Mine der Detektor- 
flache definiert ist Der Index „ i " unterscheidet also unterschiedliche Strahlenquellen- 
positionen bzw. Projektionsrichtungen. Des Weiteren ist F t der Filteroperator, der der 
Strahlenquellenposition bzw. der Projektionsrichtung der Projektion Pi zugeordnet ist, 
und pi ist die entsprechende gefilterte Projektion. 

Zum besseren Verstandnis der folgenden Schritte wird eine mathematische Beschrei- 
bung eingeschoben, die zu einer Bestimmungsgleichung fur die FUteroperatoren F t 
fiihrt. Dabei wird diese mathematische Beschreibung fur einen Computertomographen 
durchgetuhrt, ohne auf diesen beschrankt zu sein. Diese Beschreibung ist analog auf alle 
Verfahren ubertragbar, mit denen ein Projektionsdatensatz eines Untersuchungsbereichs 
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durch Akquirieren von Projektionen aus unterschiedlichen Projektionsrichtungen, wie 
es bspw. bei dem PET oder dem SPECT der Fall ist, erzeugt werden kann. Zur Ober- 
tragung auf derartige Verfahren muss im Weiteren ,,StraMenqueUenposition" durch 
,,Projektionsrichtung" ersetzt werden. AuBerdem wurde bei einem PET oder SPECT 
5 naturlich keine Verteilung von Schwachungswerten, sondern eine Verteilung von 
Aktivitatswerten rekonstruiert werden. Des Weiteren mussten, insbesondere bei PET- 
und SPECT-Verfahren, entsprechende, bekannte Vorwarts- und Ruckprojektions- 
Operatoren verwendet werden. 

10 Es gelten folgende Gleichungen: 

P t B t = id und (2) 



Pi=PiM, (3) 

15 wobei P { em Vorwartsproj ektions-Operator ist, der angewandt auf eine Verteilung jj. 
von Schwachungs- oder Absoiptionswerten im Untersuchungsbereich 13, die die 
Schwachung oder Absorption der Strahlen beschreiben, diese Werte entlang von 
Strahlen der i -ten Strahlenquellenposition vorwarts projiziert, so dass eine Projektion 
p t entsteht. Der Vomartsprojektions-Operatorbeschreibt also den Messprozess. Der 

20 Ruckprojektions-Operator B t projiziert eine Projektion entlang von Strahlen der i-ten 
Strahlenquellenposition zuruck in den Untersuchungsbereich. Wird nun eine Projektion 
entlang derselben Strahlen ruckprojiziert und vorwarte projiziert, so ergibt sich die 
urspriingliche Projektion, daher gUt Gleichung (2), wobei id ein Identitatsoperator ist 
DasheiBt, id angewandt auf eine Projektion verandert diese nicht. 

25 

Die gefilterte Ruckprojektion kann durch 



(4) 
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beschrieben werden, wobei JX das rekonstruierte Bild und N die Anzahl an Strahlen- 
queUenpositionen ist, an denen Messwerte bzw. Projektionen akquiriert worden sind, 
die fur die Rekonstruktion verwendet werden. 

Die Aufgabe der Rekonstruktion besteht darin, ein CT-Bild mit Hilfe der gemessenen 
Projektionen so zu erzeugen, dass eine Anwendung eines Vorwartsprojektions- 
Operators auf die rekonslruierten Schwachungswerte des CT-Bildes fur jede Strahlen- 
quellenposition zu den ursprunglich gemessenen Projektionen fuhrt. Ist diese 
Bedingung hinreichend erfullt, so wurde die ScMchungswert-Verteilung des Objekts 
exakt rekonstruiert. Ein Filter, der diese Bedingung zumindest im wesentUchen erfullt, 
ist ein optimaler Filter im Sinne der Erfindung. Das bedeutet, folgende Gleichung sollte 
fur jede Stranlenquellenposition j erfullt sein 



PjM = P r 



(5) 



15 



Dielndizes i undy konnenjeweils jede Steahlenquellenpositionbezeichnen, die fur die 
Rekonstruktion verwendet wird. Sie warden unterschiedlicb gewahlt, urn unten in der 
Beschreibung eine erste Stranlenquellenposition bzw. Projektionsrichtung j von einer 
zweiten Stranlenquellenposition bzw. Projektionsrichtung i unterscheiden zu konnen. 

20 

Es sei darauf hingewiesen, dass diese mathematische Beschreibung davon ausgeht, dass 
ideale Projektionen bzw. Messwerte vorliegen, d.h. Messwerte, die bspw. nicht durch 
Rauschen oder Strahlaufhartung (engl. „beam hardening") verandert worden sind Fur 
diese idealen Messwerte gilt Gleichung (5) exakt. Sind die Messwerte bspw. durch 
25 Rauschen oder Strahlaufhartung verandert, wie es bei realen Messwerten der Fall ist, so 
ist Gleichung (5) selbst fur einen optimalen Filteroperator nur im wesenllichen erfullt. 

Ein Einsetzen der Gleichung (3) in die Gleichung (5) resultiert in 
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P jjfL B < F tPi=Pj . (6) 
und Einsetzen der Gleichung (3) in die Gleichung (6) fiihrt zu 

'^ZW^-V. (7) 

Da die Gleichung (7) fur jede beliebige Verteilung // von Schwachungswerten gilt, 
mussen auch die auf diese Verteilung angewendeten Operatoren gleich sein: 
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Diese Gleichung ist aqui valent zu 




15 ScMieBlich fiihrt das Einsetzen der Gleichung (2) in Gleichung (9) zu folgender 
Besthnmungsgleichung fur die Filteroperatoren F t : 

it$Ij.t= jfllWWj = ^ . (10) 

20 Jeder Filteroperator F t , der die Gleichung (10) zumindest naherungsweise erfullt, ist 
erfindungsgemaB anwendbar. 

Im Rahmen der Erfindung wird I } _ . = P } B t F t P t B } als Iterationsoperator bezeichnet, 
damit die Bezeichnung des Operators I J t kompatibel ist nut der unten beschriebenen 
25 speziellen Avisfuhrungsform, die Filteroperatoren F t , die die Gleichung (10) erfullen, 
iterativ ermittelt. 
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Da die Ruckprojektions- und Vorwartsprojektions-Operatoren ausschlieBlich von der 
Projektionsgeometrie, d,h. von dem Verlauf der Strahlen ausgehend von der jeweiligen 
Strahlenquellenposition oder, was das Gleiche ist, von der jeweiligen Strahlenquellen- 
position, der Form des Strahlenbundels und von der Anordnung der Detektorelemente 
der Detektoreinheit, dagegen nicht von den Messwerten selbst abhangen, mussen die 
Filteroperatoren fur eine bestimmte Projektionsgeometrie nur einmal bestimmt werden. 
Bei weiteren Messungen mit der gleichen Projektionsgeometrie konnen die bereite 
ennittelten Filteroperatoren verwendet werden. 

Jedes Verfahren, das zur Ermittlung der Filteroperatoren F ( , die die Gleichung (10) 
erfullen, geeignet ist, ist erfindungsgemaB anwendbar. Zum Beispiel konnen die 
Operatoren P t B Jt PjB { und F t jeweils als Matrix dargestellt werden, wobei die 
Operatoren P t Bj und P } B % bzw. deren Matrixelemente durch emfache geometrische 
Beirachtungen aus der gegebenen Projektionsgeometrie ermittelt werden kSnnen. Bei 
15 einer derartigen Darstellung der Operatoren Uefert Gleichung (10) ein liueares 

Gleichungssystem, in dem nur die Matrixelemente des Filteroperators F { uiibekannt 
sind. Dieses Gleichungssystem kann mit bekannten Verfahren gelost werden. 

Alternativ kann Gleichung (10) iterativ mittels folgender Iterationsgleichung gelost 
20 werden: 

F ** W*W • 

Hierbei ist F* der Filteroperator der j -ten Strahlenquellenposition, der dem A: -ten 
25 Iterationsschritt zugeordnet ist. Dementsprechend ist Fj +i der Filteroperator der j -ten 
Strahlenquellenposition, der dem (*+l) -ten Iterationsschritt zugeordnet ist. Als 
Startoperator F,° fur die Filteroperatoren kann bspw. jeweils ein Identitatsoperator 
verwendet werden. 
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Im Weiteren wird das Ermitteln des fur die jeweilige Projektionsgeometrie optimalen 
Filters mit Hilfe des in Fig. 3 dargestellten Ablaufdiagramms, das sich im wesentHchen 
nach Gleichung (1 1) richtet, dargestellt Dabei wird jeder Operator als Matrix 
dargestellt. 

5 

Im Schritt 201 wird fur jeden Filteroperator Fj initial eine Matrix, bspw. eine Einheits- 
matrix, bereitgestellt. Die Matrizen fur die Ruck- und Vorwartsprojektionen konnen fur 
die jeweilige Projektionsgeometrie durch einfache geometrische Oberlegungen ermittelt 
werden. Beim Computertomographen kann im einfachsten Fall, wenn der zu rekonstru- 
10 ierende Teil des Untersuchungsbereichs (field of view - FOV) in mehrere Voxel 

eingeteilt ist und jeder Voxel einen Schwachungswert enthalt, der Vorwartsprojektions- 
Operator die Schwachungswerte derjenigen Voxel, die von einem Strahl einer Projek- 
tion durchsetzt werden, zu einem Messwert aufaddieren. Die dazugehorige Matrix 
wurde dann an den entsprechenden Stellen Einsen und ansonsten Nullen aufweisen. Der 

15 Riickprojektions-Operator kann dann bspw. so ausgebildet sein, dass ausgehend von 
einem Messwert einer Projektion aufalle Schwachungswerte, deren Voxel von dem 
Strahl, der den Messwert verursacht hat, durchsetzt worden ist, der jeweiUge Messwert 
anteilig addiert wird. Die dazugehorige Matrix besteht auch aus Einsen und Nullen, die 
so in der Matrix angeordnet sind, dass bei Anwendung der Matrix auf eine Projektion 

20 die oben beschriebene Addition von Messwerten auf die entsprechenden Schwachungs- 
werte durchgefuhrt wird. Die Erfmdung ist natiirUch nicht auf die beiden einfachen 
Beispiele eines Ruck- und Vorwartsprojektions-Operators beschrankt. 

Alternativ konnten auch Operatoren definiert werden, die eine Riick- und eine 
25 Vorwartsprojektion nacheinander ausfuhren (siehe bspw. Gleichung (19)). Diese 
Operatoren konnten anstelle eines Riick- und eines Vorwartsprojektions-Operators 
verwendet werden. Diese Operatoren bilden eine Projektion aus einem Projektionsraum 
in denselben Projektionsraum ab. 



30 
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Nachdem im Schritt 201 die Ruck- und Vorwartsprojektiora-Operatoren bestimmt und 
als startwerte Filteroperatoren F\ 0 vorgegeben worden sind, wild in den Schritten 202 
bis 207 ein Iterationsschritt ausgefuhrt. Dazu wird im Schritt 202 ein Iterationszahler k , 
der bspw. initial gleich Null ist, ran Eins erhoht 

5 

Im Schritt 203 wird eine Strahlenquellenposition j ausgewahlt, fur die in dem 
aktuellen Iterationsschritt (fc+1) nochkeinFilteroperator ermittelt worden ist 
Dies kann beim ersten Iterationsschritt (k = 1) bspw. die zeitiich erste Strahlenquellen- 
position sein. 

10 

Im Schritt 204 wird fur die Strahlenquellenposition /' und jede Strahlenquellenposition 
i jeweils ein Iterationsoperator 

1% -PjBtfPJj (12) 

15 

ermittelt, bei dem Ff der Filteroperator der i -ten Strahlenquellenposition ist, der im 
vorherigen Iterationsschritt ermittelt worden ist Wird der Schritt 204 erstmals 
durchlaufen, so ist der Filteroperator F* gleich dem initial in Schritt 201 vorgegebenen 
Filteroperator. 

20 

Die i -te Strahlenquellenposition entspricht der zweiten Projektionsrichtung und die j - 
te Strahlenquellenposition entspricht der ersten Projektionsrichtung nach Anspruch 1 . 

Im Schritt 205 wird uber alle im Schritt 204 berechneten Iterationsoperatoren gemittelt 
25 Der Mittelwert kann bspw. als arithmetischer Mittelwert bestimmt werden: 
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Alternativ konnten andere Mittelwerte verwendet werden. Aufierdem konnte, wenn 
angenommen wiirde, dass kontinuierlich Messwerte akquiriert werden, dass also keine 
diskreten SttaMenqueUenpositionen vorliegen, die Berechnung des Mittelwertes mit 
einer Summengleichung gemaB Gleichung (13) durch eine entsprechende Ihtegral- 
5 gleichung ersetzt werden. 

Der Filteroperator Ff +1 der j -ten Strahlenquellenposition des aktuellen, (k + 1) -ten 
Iterationsschrittes kann dann mittels der Gleichung 

10 Ff*=F*+id-I* (14) 

im Schritt 206 berechnet werden. Das heiBt, von der Summe aus dem entsprechenden 
Filteroperator des vorherigen Iterationsschrittes Ff, der beim ersten Durchlaufen dieses 
Schrittes gleich dem initialen Filteroperator ist, und dem Identitatsoperator wird der 
15 zuvor im Schritt 205 berechnete mittlere Iterationsoperator subtrahiert 

hn Schritt 207 wird iiberpruft, ob in diesem Iterationsschritt bereits fur jede Strahlen- 
quellenposition j ein Filteroperator Ff +1 ermittelt worden ist. Wenn dies nicht der Fall 
ist, so wird mit Schritt 203 fortgefahren. Ansonsten folgt der Schritt 208. 

20 

In Schritt 208 wird iiberpruft, ob ein Abbruchkriterium fur die Bestimmung der 
Fateroperatoren erffillt ist Das Abbruchkriterium kann bspw. jedes Mafi verwenden, 
dass die Abweichung der Filteroperatoren des aktuellen Iterationsschrittes von den 
Filteroperatoren des vorherigen Schrittes bestimmt. Ist diese Abweichung kleiner als ein 
25 festgelegter SchweUwert s , so kann das Iterationsverfahren abgebrochen werden. Das 
Abbruchkriterium kann daher bspw. folgendermafien Aussehen: 
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Hierbeisind (F/ +1 ), bzw. (F/), das t -te Matrixelement des jeweiligen Filteroperators. 
Des Weiteren konnte das Abbrachkriterium auch eine vorgebbare Anzahl an Iterations- 
schritten aufweisen, wobei die iterative Bestimmung derFilteroperatoren abgebrochen 
wird, wenn diese Anzahl erreicht oder uberscbritten ist. 

5 

Ist das Abbrachkriterium erfullt, so endet die Bestimmung der Filteroperatoren im 
Schritt 209. Ansonsten wird im Schritt 202 mit einem weiteren Iterationsschritt 
fortgefahren. 

10 Nachdem die Filteroperatoren ermittelt und die Projektionen gemaB Gleichung (1) 
gefiltert worden sind (Schritt 103 in Fig. 2), kann das erfindungsgemaBe Computer- 
tomographie-Vertahren mit Schritt 104 fortgefuhrt werden. 

Dazu wird im Schritt 104 ein Voxel F(x) innerhalb des FOV bestimmt. Dann wird im 
15 Schritt 105 eine Projektion ausgewahlt, die noch nicht zur Rekonstmktion des Voxels 
V(x) herangezogen worden ist. Verlauft kein Strahl der Projektion mittig durch den 
Voxel V(x) , so wird ermittelt, an welcher Stelle ein mittiger Strahl auf die Detektor- 
flache getrotfen ware. Der dazugehorige Messwert wird dann durch Interpolation der 
Messwerte von benachbarten Strahlen berechnet Der Messwert, der dem den Voxel 
20 mittig passierenden Strahl der Projektion zugeordnet werden kann, oder der 

entsprechende durch Interpolalion gewonnene Messwert, wird im Schritt 106 auf den 
Voxel V(x) akkumuliert. Im Schritt 107 wird gepriift, ob alle Projektionen betrachtet 
worden sind. Ist dies nicht der Fall, so verzweigt das Ablaufdiagramm zum Schritt 105. 
Ansonsten wird im Schritt 108 gepriift, ob alle Voxel V(±) im FOV durchlaufen sind. 
25 Ist dies nicht der Fall, so wird mit Schritt 1 04 fortgefahren. Wenn dagegen alle Voxel 
V(x) im FOV durchlaufen worden sind, so ist die Verteilung der Schwachungswerte, 
d.h. bspw. der Absorptionswerte, im gesamten FOV, ermittelt worden, und das 
erfindungsgemaBe Computertomographie-Verfahren istbeendet (Schritt 109). 
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Verlauft die Trajektorie, auf der sich die Strahlenquelle S relativ zum Untersuchungs- 
bereich 13 bewegt, so, dass die OperatorkombinatioB PjB t fur alle Strahlenquellen- 
positionen i und j , die auf der Trajektorie den gleichen Abstand aufweisen, gleich ist, 
so vereinfacht sich Gleichung (1 1) zu 

5 

F M =F k + Kf~IW^A • (16) 
N i 

Das heiBt, jede Strahlenquellenposition weist den gleichen Filteroperator auf, so dass 
die Schritte 203 und 207 entfallen konnen und j bspw., wie in Gleichung (16), gleich 
1 0 Eins gesetzt werden kann. Dies reduziert den Rechenaufwand fur die Bestimmung der 
Filteroperatoren erheblich. Diese mogliche Vereinfechung tritt insbesondere bei kreis- 
oder hehxformigen Trajektorien auf. 

In dem vorstehendbeschriebenen Ausfuhrungsbeispiel sind die Operatoren durch 
15 Matrizen dargestellt worden. Die Filterung einer Projektion p mittels eines Filter- 
operators F kann daher in bekannter Weise durch folgende Gleichung beschrieben 
werden: 



/=1 

Hierbei ist m die Anzahl an Messwerten pro Projektion, also bspw. die Anzahl an 
Detektorelementen, und (Fp) r ist einer der gefilterten Messwerte der Projektion p . 
Alternativ konnten die Operatoren allerdings auch analytisch dargestellt werden: 



25 



{Fp){f) = jXMpWW V^gSH 2 und (18) 

* 2 
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SR 1 

Hier ist (FpXfi) eine gefilterte Projektion an der Stelle f auf der Detektorflache, 
p(<j>) eine ungefilterte Projektion an der Stelle <f> auf der Detektorflache und fty', $) 
5 eine Filterkernel-Funktion, die in Gleichung (1 8) die Filterung beschreibt. Des Weiteren 
bezeichnen (PfijpW) eine Projektion, die zunachst riickprojiziert und danach 
vorwarts projiziert worden ist, an der Stelle </>' auf der Detektorflache und y t J W, f) 
eine Projektionskernel-Funktion. Die entsprechende Iterationsgleichung fur diese 
analyfochen Operatoren zin- Besting f(<f>',4>) lautet: 

10 

W,t) = // <*', <P) + W ~ <t>) " J K, ffl// <& 7)r,y (7, ^7 • (20) 

Hierbei ist 8ij>' - <f>) die Delta-Funktion. Die Projektionskernel-Funktion konnte im 
Schritt 201 dutch einfache geometrische "Oberlegungen aufgrund der vorgegebenen 

1 5 Projektionsgeometrie berechnet werden. Im Schritt 204 wiirden dann unterschiedliche 
Iterationsoperatoren / „ = J jrjj if, <f )// (6 riYij iV, M&V ennittelt werden, 
wonach im Schritt 205 ein Mittelwert J, durch Mittelung fiber Iterationsoperatoren mit 
vorgegebenem festen j berechnet wird. Im Schritt 206 wird dann die Filterkernel- 
Funktion if, f) des (k + 1) -ten Iterationsschrittes fur die j -te Strahlenquellen- 

20 position gemafl Gleichung (20) durch Subtrahieren des Mittelwertes Ij von der Summe 
aus der entsprechenden Filterkernel-Funkdon des A: -ten Iterationsschrittes und der 
Delta-Funktion bestimmt 

Die analytische Berechnung der Filterkernel-Funktion hat den Vorteil, dass sie im 
25 Vergleich zur Durchfuhrung der FUteroperatoren-Bestimmung mit diskreten 

Operatoren, bspw. wie obsn beschrieben mit Mairizen, keine durch die Diskretisierung 
bedingte Interpolationen benotigt und daher exakter ist. 
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Da die Messwerte einer Projektion in der Regel auf der Detektorflache an den Stellen 
der entsprechenden Detektorelemente angeordnet sind, muss die ermittelte Filterkernel- 
Funktion f{$\ fi) diskretisiert werden. Sie kann in eine Matrix umgeformt werden, so 
dass Messwerte gemaB Gleichung (17) gefiltert werden konnen. 

5 

Das erfindungsgemaSebildgebende Verfahren istnicht auf Computertomographie- 
Verfahren beschrankt Es kann bspw. auch bei einem C-Arm-System, einem PET oder 
einem SPECT angewendet werden. 

10 Bei bekannten PET- Verfahren entstehen im Untersuchungsbereich Gamma-Quanten- 
Paare durch Annihilationsprozesse (ein Elektronen und ein Positron stoBen aufeinander 
und erzeugen dabei zwei Gamma-Quanten). Da Elektron und Positron vor der 
Annihilation relativ langsam sind, bewegen sich die Gamma-Quanten eines Paares 
aufgrund des Impulserhaltungssatzes ausgehend vom Entstehungsort bekanntlich in im 

15 wesentlichen entgegengesetzte Richtungen. Urn den Untersuchungsbereich herum ist 
ein Detektorring angeordnet, der aus dem Untersuchungsbereich austretende Gamma- 
Quanten detektieren kann. Genauer gesagt werden Koinzidenzen detektiert, d.h. es wird 
die Anzahl an Koinzidenzen pro Zeitintervall und Winkelrichtung detektiert. Eine 
Koinzidenz ist dabei das im wesentlichen gleichzeitige Auftreffen zweier Gamma- 

20 Quanten an Positionen auf dem Detelctorring, die urn im wesentiichen 1 80° zueinander 
versetzt sind. Ein Detektorelement misst mit jedem der anderen Detektorelemente 
jeweils Koinzidenzen von Gamma-Quanten, die sich entlang der durch die jeweiligen 
Detektorelemente definierten Linien bewegen. Ein bestimmtes Detektorelement misst 
also zusammen mit anderen Detektorelementen Koinzidenzen in unterschiedlichen 

25 Winkelrichtungen. Alle mit einem Detektorelement pro Zeitintervall in den unterschied- 
lichen Winkelrichtungen gemessenen Koinzidenzen konnen eine Projektion bilden 
(Dies entspricht einem von einer Strahlenquelle ausgehenden Strahlenfacherbei der 
Computertomograpbie). Die Projektionsrichtung der Projektion kann dann bspw. durch 
das Detektorelement und die Mine des Detektorrings definiert sein. Die Ruckprojektion 

30 und die Vorwartsprojektion fur eine Projektionsrichtung in der jeweiligen 
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Projektionsgeometrie konnen in bekannter Weise durchgefuhrt werden, wobei die 
Projektionsgeometrie durch die jeweilige Projektionsrichtung bzw. Anordnung des 
jeweiligen Detektorelements aufdem Detekiorring und die Geometric des PETs 
definiert sein kann. Die erfindungsgemaBe Bestimmvmg von Filteroperatoren fur die 
5 jeweilige Projektionsgeometrie und die anschlieBende gefilterte Ruckprojektion fuhren 
auch bei einem PET zu einer im Vergleich zu bekannten PET-Verfahren verbesserten 
Bildqualitat 

Bei dem beschriebenen zweidimensionalen PET-Verfahren mittels eines Detektorrings 
10 kann erfindungsgemaB entsprechend Gleichung (1 6) fur jede Projektionsrichtung der 
gleiche Filteroperator verwendet werden. Da dann nur fur eine Projektionsrichtung ein 
Filteroperator zu bestimmen ist, ist der Rechenaufwand erheblich reduziert. 

Bei den SPECT-Verfahten befinden sich im Untersuchungsbereich radioaktive Isotope, 
15 die bspw. einem im Untersuchungsbereich befindlichen Patienten in die Blutbahn 
gespritzt worden sind. Die Isotope zerfaUen unter Abstrahlung von Gamma-Quanten. 
Mittels einer Gamma-Kamera mit einer Detektorflache werden aus unterschiedlichen 
Richtungen die aus dem Untersuchungsbereich austretenden Gamma-Quanten 
gemessen, d.h. es wird die Anzahl an Gamma-Quanten gemessen, die pro Zeitintervall 
20 auf der Detektorflache auftreffen. Die Messwerte, die von der Detektorflache an einer 
Gamma-Kamera-Position akquiriert worden sind, konnen als Projection bezeichnet 
werden. Die Projektionsrichtung einer Projektion kann bspw. die Richtung sein, die 
durch die Mitte des Untersuchungsbereichs und durch die Gamma-Kamera verlauft Das 
Ruck- und das Vorwartsprojizieren fur eine durch die Geometrie des SPECT und die 
25 Projektionsrichtung bestimmte Projektionsgeometrie sind bekannt, so dass das 

erfindungsgemaBe Bestimmen von Filteroperatoren und das gefilterte Ruckprojizieren 
unter Verwendung dieser Filteroperatoren auch mit einem SPECT durclrfuhrbar ist. 
Derart rekonstruierte Bilder zeigen im Vergleich zu bekannten SPECT-Verfahren eine 
signifikant bessere Bildqualitat. 

30 
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Vor der Rfickprojektion der durch PET- oder SPECT-Verfahren erzeugten Messwerte 
kann aufbekannte Weise eine Absorpbonskorrektur der Messwerte mittels 
Transmissionsmessungen durchgefubxt werden. 

5 Im Vergleich zu bekannten bildgebenden Verfahren hat die Erfindung den Vorteil, dass 
samtliche Messwerte und die exakte Projektionsgeometrie fur die Rekonstruktion 
verwendetweidenk6nnen.Sozeigtdaserfo^ 

in bekannten ComputertomograpMe-Verfahren dadurch hervorgerufen werden, dass der 
Kegelwinkel in der Rekonstruktion vernachlassigt wird. Der Kegelwinkel eines Strahles ist 

10 der Winkel, den dieser Strahl mit einer Ebene, die senkrecht zur Rotationsachse 14 
orientiert ist, einscWiefit. Ein weiterer Vorteil der Erfindung Uegt darin, dass bei einem 
Computertomographie-Verfahren mit helixformiger Trajektorie die Bildqualitat des 
rekonstruiertenCT-BildesunabhangigvondemTischvorschub pro Rotation sehr gut ist In 
bekannten Computertomographie-Verfahren ist eine Rekonstruktion von CT-Bildern mit 

15 sehr guter BildquaUtat nur fur bestimmte Tischvorschiibe pro Rotation moglich. 
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T*P7T TfrRTF.TCHENLISTE 

a Of&ungswiakel 

max 

id Identitatsoperator 

5 i 9 j Projektionsrichtungen 

jt 3 k + 1 Iterationsschritt 

p Projektion 

F i9 Fj Filteroperator 

initialer Filteroperator 

10 F * Filteroperator der i -ten Projektionsrichtung und des k -ten 
Iterationsschrittes 

Ff+ l Filteroperator der j -ten Projektionsrichtung und des (k + 1) -ten 

Iterationsschrittes 

j . Iterationsoperator 

15 /* . iterationsoperator des k -ten Iterationsschrittes 

J) Mittelwert des Iterationsoperators des k -ten Iterationsschrittes 

S Strahlenquelle 

1 Gantry 

2, 5 Motor 

20 3 KolHmatoranordnung 

4 Strahlenbiindel 

7 Steuereinheit 

I o Rekonstruktionseinheit 

II Monitor 

25 13 Untersuchungsbereich 

14 Rotationsachse 

16 Detektoreinheit 
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PATENTANSPRTJCHF. 



1. Bildgebendes Veriahren mit den Schritten: 

a) Erzeugen eines Projektionsdatensatzes eines zu rekonstruierenden Unter- 
suchungsbereichs (13) durch Akquirieren von Projektionen O y ) ausunter- 
schiedlichenProjektionsrichtimgen (J ), wobeijederProjektion O y ) eine 
Projektionsrichtung (J ) und jeder Projektionsrichtung (J ) eine Projektions- 
geometrie zugeordnet ist, 

b) Bereitstellen eines Iteiationsoperators (I jJt ) , der abhangig ist von einer ersten 
Projektionsrichtung (J ) , einer zweiten Projektionsrichtung (?) und einem der 
zweiten Projektionsrichtung (i) zugeordneten Filteroperator (F t ) , und 
Ermittelnjeweils eines Filteroperators (F t ) fur jede Projektionsrichtung (0, 
wobei die Filteroperatoren (F.) so angepasst werden, dass eine Mittelung iiber 
alle Iterationsoperatoren (F u ) , die dieselbe erste Projektionsrichtung (J ) aber 
unterscMedliche zweite Projektionsrichtungen (0 aufweisen, im wesentlichen 
einen Identitatsoperator {id) ergibt, und wobei bei Anwendung eines Iterations- 
operators (Ijj) auf eine Projektion (jpj) folgende Operationen ausgefuhrt 
werden: 

- Ruckprojizieren der Projektion (jp y ) in einer Projektionsgeometrie, die der 
ersten Projektionsrichtung (J ) zugeordnet ist, 

- Vorwartsproj izieren der ruckprojizierten Werte in einer 
Projektionsgeometrie, die der zweiten Projektionsrichtung (i) zugeordnet 
ist, 
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- Filtem der vorwarts projizierten Werte mit einem Filteroperator der zweiten 
Projektionsrichtung (i), 

- Ruckprojizieren der gefilterten Werte in einer Projektionsgeometrie, die der 
zweiten Projektionsrichtung (i) zugeordnet ist, 

5 - Vorwartsprojizieren der gefilterten, riickprojizierten Werte in einer 

Projektionsgeometrie, die der ersten Projektionsrichtung (J ) zugeordnet ist, 
c) Rekonstruieren eines Bildes des Untersuchungsbereichs aus den Projektionen 
(p.) mitHilfe einer gefilterten Ruckprojektion, wobei aufjede Projektion (pj) 
derFilteroperator (F y ) der entsprechenden Projektionsrichtung (J) 
10 angewendet wird. 

2. Bildgebendes Verfahren nach Anspruch 1, 
dadurch gekey ^^aichnet. 

dass fur jede Projektionsrichtung (j ) vor der Rekonstruktion jeweils ein Filteroperator 
15 (Fj) iterativermitteltwird, 

wobei fur jede Projektionsrichtung (J ) jedem Iterationsschritt (k+1) ein 
. . . . Filteroperator zugeordnet ist, 

wobei mehrerelterationsschritte (k+1) ausgefuhrt werden, bis ein Abbruch- 

kriterium erfullt ist, und 
20 wobei in einem Iterationsschritt (fc+1) alle Filteroperatoren die dem 

Iterationsschritt (k + 1) zugeordnet sind, ermittelt werden, bevor der nachste 
Iterationsschritt ausgefuhrt wird, und 

dass, nachdem ein Identitatsoperator (id) und fur jede Projektionsrichtung <j ) ein 
initialer Filteroperator (FJ) bereitgestellt worden sind, in einem Iterationsschritt 
25 (k + 1) zum Ermitteln des Filteroperators einer bestimmten Projektions- 

richtung (y), der dem jeweiligen Iterationsschritt (k+1) zugeordnet ist, folgende 
Schritte durchgefuhrt werden: 
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Bereitstellen eines Iterationsoperators (J*. ) , der angewandt auf eine 

Projekfcion (pj) folgende Operationen ausfuhrt: 

- Ruckprojiideren der Projektion in einer Projektionsgeometrie, die der 
bestimmten Projektionsrichtung (J ) zugeordnet ist, 

- Vorwartsprojizieren der riickproj izierten Werte in einer 
Projektionsgeometrie, die der zweiten Projektionsrichtung (i) 
zugeordnet ist, 

- Filtern der vorwarts projizierten Werte mit dem Filteroperator (F t * ) der 
zweiten Projektionsrichtung (i) , der dem vorangehenden Iterationsschritt 
(k) zugeordnet ist, 

- Riickprojizieren der gefilterten Werte in einer Projektionsgeometrie, die 
der zweiten Projektionsrichtung (i) zugeordnet ist, 

- Vorwartsprojizieren der gefilterten, riickprojizierten Werte in einer 
Projektionsgeometrie, die der bestimmten Projektionsrichtung (J ) 
zugeordnet ist 

Ermittelnjeweils eines Iterationsoperators (/*,) fiirjede Projektions- 
richtung, wobei die zweite Projektionsrichtung (?) des jeweiligen Iterations- 
operators (Ijj) die jeweilige Projektionsrichtung ist, 

— k 

Ermitteln eines mittleren Operators ( I j ) durch Mitteln iiber alle im Schritt 
ii) ermittelten Iterationsoperatoren (/*,), 

Ermitteln des Filteroperators der bestimmten Projektionsrichtung 

(J ) , der dem aktuellen Iterationsschritt (k + 1) zugeordnet ist, durch 

— k 

Subtrahieren des mittleren Operators ( Ij ) von der Summe aus dem Filter- 
operator (F*) der bestimmten Projektionsrichtung (J ), der dem voran- 
gehenden Iterationsschritts (k) zugeordnet ist, und dem Identitatsoperator 

(id). 
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3. Bildgebendes Verfahren nach Anspruch 1, 
dadurch gekennzeichnet 

dass das Erzeugen eines Projektionsdatensatzes mit einem Computertomographen im 
Schritt a) folgende Schritte aufweist: 
5 - Erzeugen eines den Untersuchungsbereich (13) durchsetzenden StraMenbundels 
(4) mit einer Strahlenquelle (S), 

- Erzeugen einer Relatrvbewegung zwischen der StrahlenqueUe (S) einerseits und 
dem Untersuchungsbereich (13) andererseits, die eine Rotation urn eine 
Rotationsachse (14) umfasst, wobei die StrahlenqueUe (S) unterschiedUche 

10 Strahlenquellenpositionen einnimmt, und 

- Akquirieren von Messwerten, die von der Ihtensitat von Strahlen in dem 
Strahlenbundel (4) jenseits des Untersuchungsbereichs (13) abhangen, mit einer 
Detektoreinheit (1 6) wahrend der Relatrvbewegung an unterschiedlichen 
Strahlenquellenpositionen, wobei Messwerte, deren Strahlen von derselben 

15 Strahlenquellenposition ausgehen, eine Projektion (j>j) bilden, jeder Projektion 

(Pj) eineProjektionsrichtung 0') zugeordnet ist, die von der jeweiligen 
Strahlenquellenposition zu einer Mitte der Detektoreinheit weist, und die 
Projektionsgeometrie, die der jeweiligen Projektionsrichtung 0' ) zugeordnet 
ist, durch die jeweilige Strahlenquellenposition und die Position der 

20 Detektoreinheit und den dadurch bestimmten Strahlenverlauf definiert ist. 

4. Bildgebendes Verfahren nach Anspruch 3, bei dem eine Trajektorie, auf der sich die 
Strahlenquelle (5) relativ zum Untersuchungsbereich (13) bewegt, so verlauft, dass ein 
Operator, der angewandt auf eine Projektion diese entlang von Strahlen einer ersten 

25 Strahlenquellenposition ruckprojiziert und die riickprojizierten Werte entlang von 
Strahlen einer zweiten Strahlenquellenposition vorwarts projiziert, gleich ist fur erste 
und zweite Strahlenquellenpositionen, die entlang der Trajektorie den gleichen Abstand 
aufweisen, und insbesondere die Trajektorie kreis- oder helixfdrmig ist, dadurch 
gekennzeichnet, dass fur jede Projektionsrichtung derselbe Filteroperator verwendet 

30 wird. 
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5. Bildgebendes Verfahren nach Anspruch 1, 
dadurch ffekerm zeichnet 

dass der Projektionsdatensatz im Schritt a) mit einem Positronen-Emissions- 
5 Tomographen oder einem Einzel-Photon-Emissions-Computertomographen erzeugt 
wird. 

6. Bildgebendes Verfahren nach Anspruch 1, 
dadurch gekennzeichnet 

1 0 dass zumindest einige der Operatoren diskretisiert sind. 

7. Bildgebendes Verfahren nach Anspruch 1, 
dadurch gekennzeichnet 

dass zumindest einige der Operatoren lineare Funktionale sind. 

15 

8. Computertomograph zur Durchfuhrung des Verfahrens nach Anspruch 3 mit 

- einer Strahlenquelle (S) zum Erzeugen eines einen Untersuchungsbereich (13) 
durchsetzenden Strahlenbundels (4), 

- einer Antriebsanordnung (2, 5) zum Erzeugen einer Relativbewegung zwischen 
20 der Strahlenquelle (S) einerseits und dem Untersuchungsbereich (13) 

andererseits, die eine Rotation urn eine Rotationsachse (14) umfasst, 

- einer Detektoreinheit (1 6) zum Akquirieren von Messwerten, die von der 
Intensitat von Strahlen in dem Strahlenbiindel (4) jenseits des 
Untersuchungsbereichs (13) abhangen, wahrend der Relativbewegung, 

25 - einer Rekonstruktionseinheit (10) zum Rekonstruieren eines 

Computertomographie-Bildes des Untersuchungsbereich (13) aus den 
Messwerten, 

- einer Steuereinheit (7) zum Steuem der Strahlenquelle (S), der 
Antriebsanordnung (2, 5), der Detektoreinheit (16) und der 

30 Rekonstruktionseinheit (10) entsprechend den Schritten nach Anspruch 3. 
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9. Computerprograrnm fur eine Steuereinheit (7) zur Steuerung einer Strahlenquelle (S), 
einer Detektoreinheit (1 6), einer Antriebsanordnung (2, 5) und einer Rekonstruktions- 
einheit (10) eines Computertomographen zur Durchrohrung der Schritte nach Anspruch 
5 3. 
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ZUSAMMENFASSUNG 

Bildgebendes Verfahren mit gefilterter Riickprojektion 

Die Erfindung betrifft ein bildgebendes Verfahren, bei dem ein Projektionsdatensatz 
eines zu rekonstruierenden Untersuchungsbereichs durch Akquirieren von Projektionen 
5 aus unterschiedlichen Projektionsrichtungen, insbesondere mit Hilfe eines Computer- 
tomographen, erzeugt wird. Fur jede Projektionsrichtung wird ein Filteroperator 
ermittelt, der an eine der jeweiligen Projektionsrichtung zugeordneten Projektionsgeo- 
metrie optimal angepasst ist. Aus dem Projektionsdatensatz wird ein Bild des Unter- 
suchungsbereichs durch Filtern der Projektionen mit den ermittelten Filteroperatoren 
1 0 und durch Riickprojizieren dieser gefilterten Projektionen rekonstruiert. 



Fig. 3 
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